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The dynamic electric dipole polarizability function for the magnesium atom is formed by assem¬ 
bling the atomic electric dipole oscillator strength distribution from combinations of theoretical and 
experimental data for resonance oscillator strengths and for photoionization cross sections of valence 
and inner shell electrons. Consistency with the oscillator strength (Thomas-Reiche-Kuhn) sum rule 
requires the adopted principal resonance line oscillator strength to be several percent lower than 
the values given in two critical tabulations, though the value adopted is consistent with a number 
of theoretical determinations. The static polarizability is evaluated. Comparing the resulting dy¬ 
namic polarizability as a function of photon energy with more elaborate calculations reveals the 
contributions of inner shell electron excitations. The present results are applied to calculate the 
long-range interactions between two and three magnesium atoms and the interaction between a 
magnesium atom and a perfectly conducting metallic plate. Extensive comparisons of prior results 
for the principal resonance line oscillator strength, for the static polarizability, and for the van der 
Waals coefficient are given in an Appendix. 

PACS numbers: 34.20.Cf, 32.10.Dk, 31.10.-l-z 


I. INTRODUCTION 

Magnesium is an abundant element currently of interest in several applications. Analysis of photo association 
spectroscopy for the Mg dimer mm indicates that the s-wave scattering length for collisions between two ground 
state Mg atoms is positive [1], with the accuracy of the determination affected by the remnant uncertainty in the 
value of the atom-atom van der Waals constant [5] . (A Bose-Einstein condensate of Mg atoms has not been created 
experimentally, to date.) However, a determination of the leading term in the long-range interaction of electronically 
excited states of the Mg dimer, which is related to the principal resonance line oscillator strength, using molecular 
spectroscopy has been elusive [3]. Mg atoms were investigated as a possible sympathetic cooling agent in collisions 
with NH [U |S] and Lonij et al. [5] theoretically explored the interaction of an Mg atom with a wall for applications to 
atom interferometry. And, while Mg is abundant in the solar system, it is interesting to note that the Mg principal 
resonance line was recently detected in the observation of the exoplanet HD 209458b using transit spectroscopy [7]. 
Modeling Mg absorption in exoplanet atmospheres depends proportionally on the principal resonance line absorption 
oscillator strength [5] , the value of which in turn affects the use of that line as a probe of escaping atoms in exoplanet 
atmospheric spectroscopy [9]. 

An analysis of the dynamic electric dipole polarizability function of Mg is valuable for several reasons. First, because 
the principal resonance line oscillator strength is an important contributor to the function, it is possible to determine 
a value that achieves consistency with oscillator strength sum rules. Secondly, it is desirable to have an independent 
assessment of the completeness of existing elaborate calculations of the function itself, for which extensive tabulations 
from two different types of calculations are available [iniiii]. 
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Calculations of dynamic electric dipole polarizabilities are of intrinsic theoretical interest due to the challenges 
inherent in treating correlations and excitations of all electrons quantum-mechanically at different photon energies m 
Such calculations are necessarily important benchmarks for theoretical methods applied to photoabsorption 
photodetachment m, blackbody radiation shifts [TMTS] and AC Stark shifts [H] , magic wavelengths m, and parity 
non conservation amplitudes m, as well as being helpful in the ongoing development of density functional theory 
(DFT) methods for dispersion forces (cf. P^26) b In addition, for metals experimental data at a wide spectrum of 
photon energies are relatively scarce though X-ray data, and sometimes optical data HZlllHl, are available. There 
is recent progress for systematically measuring static polarizabilities [23] ■ Many theoretical approaches are available, 
but their reliability in calculating dynamic polarizabilities can be difficult to gauge without critical evaluation, but 
critical evaluations are limited to the static polarizabilities [SOl - l^ . Nevertheless, dynamic polarizabilities are of great 
utility in calculating coefficients appearing in certain potential energies, particularly van der Waals constants (for 
investigations of ultra-cold collisions, for photo association spectroscopy, and for ultra-cold gas studies) and Lennard- 
Jones constants (for atom-surface interactions, where recent applications include tests for gravity-related new physics 
at submillimeter distances [83U35] , optical clocks [Ml [37] , atom-graphene interactions [M] , noncontact van der Waals 
friction |M|, and interactions between nanostructures [iP]'). 

In a previous paper treating the sodium atom [41j , a semi-empirical theory utilizing oscillator strength sum rules and 
input data from experiments and calculations predicted a value of the van der Waals coefficient m. which was found 
to be in harmony with subsequent experimentally determined fits from photo-association spectroscopy data [421 143j 
and ab initio theoretical methods [44] . In the case of sodium, the availability of precise measurements of the principal 
resonance line oscillator strength from photo-association spectroscopy and of the static electric dipole polarizability 
from atom interferometry augmented the semi-empirical analysis [41] . In the case of magnesium, such data are not 
available. Therefore, in this paper, for Mg, the electric dipole oscillator strength distribution is composed using extant 
data on electric dipole oscillator strengths, photoabsorption cross sections, and energies obtained experimentally and 
theoretically. As will be shown, consistency with the oscillator strength (Thomas-Reiche-Kuhn) sum rule requires 
a value of the principal resonance line oscillator strength that is several percent lower than values listed in critical 
tabulations by Morton [45] and by Kelleher and Podobedova [46]. Other evidence for the value adopted is given. The 
static electric dipole polarizability is evaluated and compared with other values. The dynamic polarizability function is 
calculated and compared with previous results obtained by Porsev et al. [10] using configuration interaction and many- 
body perturbation theory with core contributions (CI-MBPT) [17] and by Jiang at al. E] using the configuration 
interaction with semi-empirical core-valence interaction (CICP) method. The present dynamic polarizability function 
is used to evaluate the van der Waals constant, Axilrod-Teller-Muto constant, and atom-surface interactions. Results 
from the literature for the principal resonance line oscillator strength, static electric dipole polarizability, and van der 
Waals constant are collected and compared in the Appendix. 


II. DIPOLE OSCILLATOR STRENGTH SUM RULES 


The absorption oscillator strength from the ground state |0) with eigenvalue Eq to an excited state |n) with 
eigenvalue En is 
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where is the position vector of electron and N is the number of electrons. Atomic units are used throughout 
unless otherwise specified. 

Denoting by S'„ the sum-integral (the sum over all discrete transitions excluding the initial state and the integration 
over all continuum states), the resultant sum rules are 

5(0) = S'jr. = N, (2) 


with N = 12 for Mg, 


and 


5(-l) 


S'jn/iE, 



(3) 


5(-2) = 5'Jn/{En - E^f = a(0), 


(4) 


where a(0) is the static electric dipole polarizability. The dynamic electric dipole polarizability function is 

= S' (5) 

[En - EqY - 

where lo is the photon energy. By direct integration the 5(—1) sum rule is related to the atom-wall interaction 
coefficient m 


the van der Waals coefficient is 


1 

Cs = |5(-1) = — y dwa{iuj), 

Ce = — / duj [a{iuj)]'^, 

Jo 


and the Axilrod-Teller-Muto coefficient is 
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Cg = — doj [a{iuj)]^. 
^ Jo 


( 6 ) 


(7) 


( 8 ) 


III. OSCILLATOR STRENGTH DISTRIBUTION 

A magnesium atom has twelve electrons. Their configuration is (ls^2s^2p®3s^) ^5o. 


A. Discrete transitions 

Sources with tabulations of values, experimental and theoretical, for the absorption oscillator strengths are 
Mitchell |3H]; Mendoza and Zeippen [1^, Ray and Mukherjee [50], Jonsson and Fischer m. Hamonou and Hib- 
bert [52], and Derevianko and Porsev [53] . 

There are numerous theoretical determinations of the oscillator strength for the principal resonance line (3s^) — 

(3s3p) A detailed survey is given in the Appendix. Reliable theoretical calculations range from 1.709 to 1.76, 
and there were at least ten experimental determinations as of 2003 [JS]. In a critical review, Morton [35] adopted a 
value of 1.83 ± 0.03 based on a weighted mean of the ten experimental values. A long-standing discrepancy between 
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experimental and theoretical trends was noted previously [MHST]- The tendency of high-level theoretical results to be 
less than 1.8 was noted recently by Zatsariny et al. m who calculated a nonrelativistic value of 1.738 and pointed 
out that a “very extensive and essentially converged multiconfiguration Hartree-Fock (MCHF)” ab initio calculation 
by Jonsson, Fischer, and Godefroid [S5] found 1.717. Jonsson, Fishcher, and Godefroid [S5], using the observed 
transition energy to evaluate the oscillator strength, obtained 1.710 and also pointed out (see their Table 11) that 
(as of 1999) theoretical values were consistently smaller than the experimental ones (to date) by about 5 percent. 
Recently, Derevianko and Porsev [S3] quote for the matrix element governing the line strength a value 4.03 ±0.02 with 
error of 0.5% based on their calculations from 2001 [35]; the corresponding oscillator strength using the experimental 
transition energy [ISj is 1.73±0.02. The original NIST (NBS) tabulation of 1969 [33] adopts 1.8±0.18 from an average 
of the experiments of Refs. [60l[6T] and the calculation of Weiss [62]. The value 1.8 (cited as a private communication 
from A. W. Weiss) is given with 3 percent error (±0.05) at the 90 percent confidence level in the 2008 NIST revised 
tabulation [46] . 

1 adopt the value 1.75, which is the lower limit of the value 1.80(5) from Kelleher and Podobedova [33] and the 
upper limit of the value 1.73(2) recommended by Derevianko and Porsev. 

Mitchell [3H] used the anomalous dispersion (hook) method to measure the second (3s-4p) through sixth (3s-8p) 
resonance transition oscillator strengths and found, respectively, 0.107±0.0019, (2.27±0.12) x 10“^, (8.53±0.46) x 10“^, 
(4.11±0.36) X 10“^, and (2.34±0.15) x 10“^, all determined relative to the value of 1.72 for the principal transition that 
Victor and Laughlin [63] calculated using a semi-empirical model potential method. Other theoretical determinations 
of / values for the second and higher resonance transitions are those of Saraph [53] , Mendoza and Zeippen [33] , Ghang 
and Tang [65] and Zatsarinny et al. m- 1 adopt the values of Ghang and Tang, who calculated for the principal 
to sixth resonance lines, respectively, oscillator strengths 1.75, 0.111, 0.024, 0.0091, 0.0043, and 0.0024. These five 
discrete (second to sixth resonance) transitions contribute 0.151 to the S'(O) sum and 2.81 to the S{—2) sum for the 
3s shell. 

Including the principal resonance line, the discrete transition contribution from the 3s shell to S'(O) is 1.90 and to 
5(-2) is 71.4. 


B. Continuum transitions 

A number of sources exist for the continuum oscillator strengths corresponding to the ejection of a 3s electron [551 - 
155] . The photoionization cross sections calculated using a variational MGHF method by Fischer and Saha m are in 
good agreement with the experimental results of Wehlitz et al. [82]. The threshold cross section of 2.5 Mb calculated 
by Fischer and Saha m is slightly larger than both the value 2.1 ±0.3 Mb measured by Fung et al. [53] and the value 
2.36±0.02 Mb that Parkinson, Reeves, and Tomkins [53] found by extrapolation from the measured discrete oscillator 
strengths of Ref. [35]. Wehlitz et al. [53] normalized their own measurements to the value 2.1 Mb at threshold. At 
a photon energy of 30 eV, a recent calculation by Pindzola et al. [53] using a time-dependent close-coupling method 
with an effective core potential gives (T3s(30eV) = 0.217 Mb, while Verner et al. [75] calculate 0.255 Mb. At 80 eV, 
the “complete” experiment of Haussman et al. [83], for which the total absorption cross section of Ref. [87] was 
used for normalization, yields cr3s(80eV) = 0.080 ± 0.011 Mb from the main transition compared to the value 0.087 
from Verner et al. Haussman et al. also measured an additional 0.014 ± 0.004 Mb contribution from satellites. The 
measurements of Wehlitz et al. [53] are adopted for energies from threshold to 11.6 eV and the 3s cross section data 
were extended to higher energies using the results of Verner et al. m- 
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The contributions to the three sum rules S{0), S'(—1), and S{—2) from the continuum are, respectively, 0.261, 0.162, 
and 0.277. Combining these with the discrete contributions from the Sec. EnH the total valence shell contributions 
are 2.16, 11.8, and 71.7, which are listed in Table |l] The calculated valence shell contribution of 2.16 to the 5'(0) 
sum confirms that the contribution from the 3s electron to S'(O) is greater than 2 [551 IBS], indicating configuration 
interaction of the valence electrons with the core electrons. Maeder and Kutzelnig m obtained 2.06 using a model 
potential including core-valence correlation. 

The excitation and ejection of K shell electrons was considered by Verner et al. |75|, Kutzner et al. EH, and 
Hasoglu et al. [52]. Hasoglu et al. used i?-matrix methods to calculate excitation of K shell electrons to np states 
resonant below the threshold. The resonances are estimated to contribute 0.03 to the value of 5'(0). The relativistic 
random-phase approximation modified to include relaxation effects (RRPAR) calculations of Kutzner et al. are in 
good agreement with the cross sections calculated by Verner et al. from threshold to high energies. Banna et al. [55] 
measured the shake-up peak just above the threshold, but absolute cross sections are not available. The cross sections 
of Verner are adopted from the threshold for ejection of a Is electron at 1310.9 eV. The K shell contribution yields 
1.56 to S'(O) and it is negligible for S'(—1) and S{—2). The Is contribution to S'(O) found here for Mg is comparable 
to that found for Na m- 

Deshmukh and Manson m, Nasreen, Manson, and Deshmukh [51], Kutzner et al. [5T], and Verner et al. Hg 
calculated the partial cross sections for the ejection of a 2s electron. The cross section value 0.3 Mb at the threshold 
energy of 94.0 eV is adopted [511IM] and linearly joined to the results of Verner et al. at 270 eV, which are used for 
higher energies. For the 2s shell the contribution to the sums 5(0), 5(—1), and S{—2) are, respectively, 1.02, 0.1, and 
0 . 02 . 

The remaining oscillator strength must come from the 2p shell and the expected contribution to 5(0) is 7.26. 

Deshmukh and Manson HQ], Altun Nasreen, Manson, and Deshmukh [94], Kutzner, Maycock, and Thorarin- 
son m, investigated photoionization of a 2p electron. There a number of resonances corresponding to excitation to 
autoionizing states m and they contribute significantly to the oscillator strength. Measurements of the cross section 
by Haussman et al. [55] at 80 eV found that the resonances constitute 25% of the 2p shell photoionization cross 
section. The cross sections averaged over the resonances from the correlated length gauge many-body perturbation 
theory (MBPT) calculations of Altun (Figure 9 of Ref. [77]) are adopted from 63.29 eV to 344.89 eV, giving <T2p(80 eV) 
= 6.2 Mb which is slightly larger than the reference value used by Haussman et al. The calculations of Verner et al. 
are used from the threshold energy of 54.9 eV up to 63.29 eV and for energies above 344.89 eV. These data yield 
for the 2p shell the contributions to the sums 5(0), 5(—1), and 5(—2), respectively, 7.33. 1.61, and 0.45. If the 
calculations of Altun are multiplied by a factor of 0.985, the cross section at 80 eV becomes 6.1 Mb and the sums 
5(0), 5(—1), and 5(—2) are calculated to be, respectively, 7.26, 1.60, and 0.44 and the values are listed in Table|T] 

The contributions to the sums 5(0), 5(—1), and 5(—2) from the excitation of the Is, 2s, 2p, and 3s electrons are 
summarized in Table H] 


C. Discussion 

Stwalley [55] , Pal’chikov and Ovsiannikov [5^ , Ovsiannikov et al. [55] , and Sarkisov et al. [55] constructed oscillator 
strength distributions of Mg for calculations of dynamic polarizabilities, considering valence transitions. Sarkisov et 
al. [28] included an estimate of 2p excitations. 

From their tabulated data, the results of Pal’chikov and Ovsiannikov [95] and Ovsiannikov et al. [5^ indicate a total 
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TABLE I: Contributions to the sum rule S{k) for Mg. 



5 ( 0 ) 

5(-l) 

5 (- 2 ) 

Is 

1.56 

0.02 


2s 

1.02 

0.10 

0.02 

2p 

7.26 

1.60 

0.44 

3s 

2.16 

11.8 

71.7 

Total 

12.0 

13.5 

72.2 


discrete contribution of 1.9 to the <S'(0) sum rule for 3s discrete transitions in agreement with the present result. They 
used a value of 1.73 for the principal resonance line oscillator strength, which offsets in the sum rule their slightly 
larger value of 0.122 for the second resonance line oscillator strength, compared to the present adopted values of, 
respectively, 1.75 and 0.111. For the S{—2) sum, they find for the second to sixth resonance transitions a contribution 
of 3.12 to S{—2), to be compared to the present value of 2.81. The difference between their value and the present 
value is due primarily to the different values of the second resonance line oscillator strength. Their total value (valence 
electrons) for S{—2) is 71.39, while the present value for the 3s shell is 71.7. 

Including all shells, the present value of a(0) is 72.2 from the <S'(— 2 ) sum rule, see Table [I] It lies only 0.2 above 
the range of values 71.3(7) recommended by Porsev and Derevianko [47] and it is compared with a number of other 
theoretical calculations in the Appendix. 

Note that if the value 1.8 is adopted for the principal oscillator strength [35], without any other adjustments to the 
adopted data, the present 5'(0) sum becomes 12.05 and a(0) becomes 73.4, which is far beyond the value recommended 
by Porsev and Derevianko. The value 1.83 for the oscillator strength adopted by Morton |3^ is more difficult to 
reconcile within the present analysis. The S'(O) sum becomes 12.08 and the value of a(0) becomes 74.6. Sarkisov et 
al. [55] use the oscillator strength data from Morton |3S| and estimate the 2p and 3s continuum contributions using 
the data from Verner et al. ESI and find a(0) = 73.6. Stwalley’s early calculation jSS] used the value 1.82(5) for the 
principal resonance line oscillator strength (see his Ref. 7 for sources) and obtained an estimate q;(0) = 75.0 ± 3.0. 

A recent experiment |29j using a pulsed cryogenic molecular beam electric deflection method obtained a value 
a(0) = 59(15), which is not sufficiently accurate to discriminate between theoretical calculations. 

Reshetnikov et al. |89j explored the relationship between the uncertainty in a(0) and the uncertainty in the lifetime 
of the first resonance transition in two-valence electron atoms and ions. Their formalism allows a valence shell 
contribution to N that is not exactly 2, as found here and in Ref. [SDj. In terms of the the valence contribution, N^, 
the principal resonance line oscillator strength and the excitation energies of the first and second resonance 

transitions, respectively, i? 3 s, 3 p and i? 3 s, 4 p, they give 


a(0) 


/3s,3p . 

^is,3p ‘^^3s,Ap 


(9) 


and an uncertainty estimate for the polarizability 

Aa(0) 


Ne — f3s,3p 
2^3s,4p 


( 10 ) 


Using the present adopted value f 3 s, 3 p = 1-75, calculated value = 2.16, and transition energies [Mj E 3 s, 3 p = 
0.159 705 and E^s.Ap = 0.224 840, yields an estimate a(0) = 72.7 ± 4. Likewise, using their formula for estimating 
the uncertainty of / 3 s, 3 p, given the present calculated value a(0) = 72.2, yields / 3 s, 3 p = 1-79 ± 0.07. The formulae of 
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FIG. 1: The dynamic dipole polarizability function a{iu}) from the present calculations. 

Reshetnikov et al. [89| demonstrate that the present results are mutually consistent, but the estimates obtained are 
not sufficiently precise to allow selection of a particular value of /ss.sp from the many available values, see Appendix. 

The availability of a more accurate measurement of q;( 0) and a definitive measurement of the principal resonance 
line lifetime would significantly improve the present model m- Nevertheless, the values adopted here, in particular 
f 3 s, 3 p = 1.75, generate sum rules that are consistent and not in contradiction with other major studies, while a value 
of f 3 s, 3 p > 1-8 is inconsistent. 


IV. DYNAMIC ELECTRIC DIPOLE POLARIZABILITY FUNCTION 


The dynamic electric dipole polarizability function at imaginary frequencies is constructed using the discrete and 


continuum oscillator strength data as assembled in Sec. Ill 


The continuum oscillator strength distribution is given in terms of the photoionization cross section a{E) by 

(t{E) 


dE 


E > 0.281, 


27 r 2 Offs’ 

with Ofs the fine structure constant, and the dynamic dipole polarizability at imaginary energy is 




fn 


(A„ - Ao)2 + 


+ / dE 


df/dE 


A2 


• UJ 


2 ’ 


( 11 ) 


( 12 ) 


The function resulting from the analysis in Sec. Ill is shown at low energies in Fig. It may be compared 

with the calculations of Derevianko et al. [10] and those of Jiang et al [II]. The present function a{iuj) was evaluated 
at the fifty energies corresponding to the energies ujk of a 50-point quadrature, as fisted in Table A of Ref. cni, and 
the energies for a 40-point quadrature listed in Table C of Jiang et al. In Fig. [^the data are plotted. Agreement 
is very good between the present results and the CTMBPT results of Ref. m- There are noticeable discrepancies 
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photon frequency (units of e^/ag) 

FIG. 2: For the dynamic dipole polarizability function a{iu}), comparison between the present values (line), the configuration 
interaction with semi-empirical core-valence interaction (CICP) values from Ref. [TT] (circles), and the configuration interaction 
and many-body perturbation theory with core interactions (CI-MBPT) values from Ref. [10] (squares). 

between the present results and the CICP results of Ref. m- To further investigate the discrepancies, in Fig. the 
percentage difference between the values from the functions given in Ref. [TU] or Ref. [TT] and the present values 
are shown. The present model and the calculations of Ref. m agree within several percent at all energies. The 
present values are larger than those of Ref. uni at low energies, due to the larger principal oscillator strength adopted 
herein. The several percent agreement across all energies is satisfactory and a more detailed analysis might await an 
accurate experimental value for a(0) and a definitive measurement of the principal resonance line oscillator strength. 
In contrast, the differences between the present model and the model of Ref. m are made apparent in Fig. [^showing 
that the CICP model of Ref. ffH yields larger values for a{iuj) in the energy range of 1 to 200 au. The difference 
may arise due to the choice of “effective” core oscillator strengths in the CICP model [iniT]. The percent difference 
peaks at about 16 % around 5-6 e^/ag (135-160 eV) placing the missing oscillator strength of the CICP model in the 
inner s shells, where the “effective” oscillator strengths are placed to model inner shell absorption [IHIIT]. 


V. APPLICATIONS 


A. van der Waals coefficient Ce 


The long-range potential energy between two Mg atoms separated by a distance R is —Cq/R^, where Cg is given 
by Eq. 0. For the van der Waals coefficient, I find Cg = 642.4 by evaluating Eq. 0 using the quadrature method 
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photon frequency (units of e^/ao) 

FIG. 3: For the dynamic dipole polarizability function a(itu), percent difference comparison between the present values and 
the CI-MBPT values from Ref. |10| and between the CICP values from Ref. [TT| . The quantity plotted is 100 x [(other) — 
(present)]/present, where “other” is either Ref. [10] or m- The filled squares are the percent differences of the present values 
from the CI-MBPT values and the filled circles are the percent differences of the present values from the CICP values. The 
plotted values show that the present a{iuj) is within several percent of the CI-MBPT values across the range of energies and 
show that the CICP values are larger for photon energies between roughly 1 and 200 atomic units. 

of Ref. Uni, 

3 

(13) 

TT ^' 

1 = 1 

using the present values of a(iuj) evaluated at the energies LOj and with the weights Wj given in Table A of Ref. m- 
Porsev and Derevianko |98| quote accuracy of 2% or better for their value Cq = 627(12) obtained using a semi- 
empirical hybrid relativistic many-body perturbation theory (CI-MBPT) approach. The present value lies just above 
their range, mainly corresponding to their principal resonance line oscillator strength of 1.73, compared to the present 
adopted value of 1.75, thus, 642 x (1.73/1.75)^ « 627. The present result improves upon the earlier empirical estimate 
of 683(35) by Stwalley | 88 |. A more detailed survey and comparison of other determinations of Cq is given in the 
Appendix. 


B. Atom-wall coefficient C 3 

The long-range potential energy of an Mg atom at distance z from a perfectly conducting wall is —Ca/z^, where C 3 
is given by Eq. (§. Mitroy and Bromley m calculated C 3 = 1.704 using the CICP approach, while the CI-MBPT 
value is 1.666 [lOj . Lonij et al. [^ gave an approximate value of 1.51 using a limited 4-parameter model for the 
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dynamic polarizability. The value of C 3 is known to be sensitive to the completeness of the description of the core 
electrons PIIMIITMI. 

The present value is C 3 = 1.69 using C 3 = |S'(—1) and the value of <S'(—1) from Tableand 1.687 using Eq. ^ 
and the quadrature from Ref. nni, 


C3 




(14) 


Both of the present values (sum rule and quadrature) are in larger than that of Ref. [TU] . 
of C 3 from the CICP calculations, Ref. m, is consistent with the relatively larger values 
Sec. IV and shown in Figs. and 


The slightly larger value 
of a{iuj), as discussed in 


C. Other properties 


The Axilrod-Teller-Muto coefficient Cg, Eq. (|^, characterizes the mutual long-range interaction potential of three 
atoms. The value obtained by Mitroy and Bromley [57] is 33 380 and that obtained by Porsev et al. [10] is 33 241. 
Using the dynamic polarizability function I evaluated, Eq. using the quadrature of Ref. [TO] . 


Cg 


3 

TT 


(15) 


and obtained Cg = 34480. 

The larger value for Cg found here mainly reflects the larger principal oscillator strength 1.75 adopted compared 
to the principal oscillator strengths found in Refs. |97| and [l0|. The oscillator strength appears as a cubic power in 


Eq. (15) through a{iuj), see Eq. (12). For example, comparing to Ref. uni, which used a principal oscillator strength 
of 1.73, scaling the present value I obtain 34 450 x (1.73/1.75)^ = 33 310, which is within 0.2 % of the value of Ref. [10] . 


VI. CONCLUSION 

Experimental and theoretical data were assembled and used to formulate the dynamic polarizability function for 
Mg. I find that consistency in the sum rules can be achieved using the adopted value of the principal resonance line 
oscillator strength to be 1.75; lower than the curated values of 1.83 [IS] and 1.8 [15], but in agreement with theoretical 
calculations. Comparisons of the dynamic dipole polarizability functions from the present work and those calculated 
using the CI-MBPT approach and the CICP approach were presented. Good agreement (within several percent) was 
found with the CI-MBPT results over all photon energies providing an independent confirmation of the CI-MPBT 
approach for Mg m- For the CICP method the differences were more pronounced, approaching 16% at energies 
around 5-6 e^fag, or about 135-160 eV, indicating that the “effective” oscillator strengths of Refs. [97] and [TT] may 
not completely model oscillator strengths corresponding to the core electrons. To improve the present model it would 
be valuable to have more accurate experimental measurements of the polarizability and a definitive measurement 
of the principal oscillator strength. Where sufficient and reliable data is available, the present methodology can be 
applied to other atoms. Note added in proof: A recent calculation [Y. Singh, B. K. Sahoo, and B. P. Das, Phys. 
Rev. A 88, 062504 (^013)], using a relativistic coupled cluster theory with all singly and doubly excited configurations, 
found a polarizability value of 72.54(50). I thank Dr. Singh for communicating this result. 
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Appendix: Values from the literature 


In this Appendix, values for Mg of the principal resonance transition oscillator strength, of the static electric 
dipole polarizability, and of the van der Waals constant are collected from the literature. Some earlier collections 
include Refs. [iSHSOl 175] for the principal oscillator strength, Refs. [311 EOl HH HOI 1271 UOlj for the polarizability, and 
Ref. [97] for the van der Waals constant. 

For the oscillator strength, as discussed in the Sec. Ill it was noted in several recent papers [SUIT] that in general 


the most sophisticated theoretical calculations lie several percent below the published experimental values, see also 
earlier similar comments in Refs. [53], [5S], and [H]. Also, it was noted that the “best” calculations lie below m 
the adopted value of 1.8 in the NIST tabulation [46]. As shown in Table [ll| the configuration interaction (Cl), multi¬ 
configuration Dirac-Fock (MCDF), and multi-configuration Hartree-Fock with Breit-Pauli interactions (MCHF-I-BP) 
calculations, are all in general agreement; the latter two methods include relativistic effects. In addition, the semi- 
empirical model potential (MP) calculation of Victor and Laughlin [53] agrees with the MCHF calculation. In 
the Table, where treatment of core-valence correlation was included the suffix CV is appended. In addition, the 
relativistic configuration interaction with the Breit interaction (RCI-|-Breit ) of Ref. |103j are in close agreement 
with the configuration interaction Dirac-Fock with core polarization (GIDF-GP) calculation of Ref. |104j . CI-I-MBPT 
methods use complete relativistic Cl calculations for the valence electrons in a frozen core combined with MBPT to 
account for core-valence interactions. The Cl frozen core calculations of Saraph [64] and Chang and Tang [65] and 
the CI-CV calculations of Moccia and Spizzo [75] are in close agreement with values between 1.75 and 1.76. There is 
substantial theoretical evidence for a value of the principal oscillator strength around 1.75. 

Many calculations of the static polarizability Q!(0) are available. There are several good tables containing summaries 
of other earlier works. In particular, from 1976, Reinsch and Meyer m and from 1991, Archibong and Thakkar |108j . 
see also Schwerdtfeger [35] . Thakkar and Lupinetti m recommend a theoretical value of 71.22 ±0.36, which includes 
a relativistic correction of —0.35, and Chu and Dalgarno |109j recommend 71. 


Table III lists some of the results from the literature. The quantum defect theory (QDT) value from Chernov 
et al. |110j models the response of the valence electrons only. The pseudopotential (PP) [90] and model potential 
(MP) calculations |11111112j model the response of the valence electrons with inclusion of effective potentials to 
treat the core electrons. An effective core potential is used similarly in the configuration interaction core potential 
(CICP) calculation of Muller, Flesch, and Meyer |113j . The CICP calculation of Mitroy and Bromley [57] utilizes a 
model that treats core excitation using effective oscillator strengths designed to reproduce the core polarizability. The 
multi-reference configuration interaction (MRGI) calculation of Partridge et al. |114j gives 71.2. They also calculated 


a CI-GV value of 70.3 (not listed in Table III), which is in good agreement with the similar calculation of 70.9 by 
Hamanou and Hibbert [^, but the average value of 74.37 from the CI-GV calculations by Moccia and Spizzo |115j is 
significantly larger. The coupled cluster double-excitation with contributions of single and triple excitations [CCD±ST 
(CCD)] model of Castro and Canuto [116] yields a value of 70.89, somewhat lower than the MRGI value and the fourth 
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TABLE II: A comparison of results for the principal resonance transition absorption oscillator strength for Mg. Abbreviations 
for methods are defined in the text. Where a calculation was carried out in the length gauge (LG) and the velocity gauge (VG), 
the simple average is listed and the LG and VG values are listed as a table footnote. 


Method 

/ 

Source 

Ref. (Year) 

MCDF-CV 

1.709 

Jonsson and Fischer 

[SI] (1997) 

MCHF 

1.717 

Jonsson, Fischer, & Godefroid 

jSD (1999) 

MP 

1.717 

Victor and Laughlin 

[63] (1973) 

MCHF-hBP 

1.719 

Fischer, Tachiev, & Irimia 

[T02] (2006) 

RCI-|-Breit 

1.722 

Cheng et al. 

[T03| (2011) 

CIDF-CP 

1.72 

Stanek, Glowacki, & Migdalek 

[T04| (1996) 

CI+MBPT“ 

1.724 

Savukov and Johnson 

flMI (2002) 

cU 

1.725 

Mengali and Moccia 

fT06| (1996) 

CF 

1.73 

Nesbet and Jones 

[Tot] (1977) 

CI-tMBPT‘* 

1.73(2) 

Derevianko and Porsev 

[S3] (2011) 

CICP 

1.732 

Mitroy and Bromley 

[97] (2003) 

Cl 

1.735 

Hamonou and Hibbert 

[S2] (2008) 

cr 

1.737 

Weiss 

[SI] (1967) 

MCHF-CV 

1.738 

Zatsarinny et al. 

[St] (2009) 

MCHF^ 

1.747 

Fischer 

[S3] (1975) 

Cl-frozen core 

1.75 

Chang and Tang 

[SS] (1990) 

CI-CV» 

1.755 

Moccia and Spizzo 

[76] (1988) 

Cl-frozen core 

1.76 

Saraph 

[S3] (1976) 

NIST adopted'* 

1.8 

Kelleher and Podobedova 

[IS] (2008) 

Experiment’ 

1.83(3) 

Morton 

[3S] (2003) 


“Mean of LG and VG line strength with experimental transition energy from |46| 

^Mean of LG 1.72 and VG 1.73 
“Mean of LG 1.746 and VG 1.717 

^Line strength with experimental transition energy from |46| 

“Mean of LG 1.773 and VG 1.701 
^Mean of LG 1.757 and VG 1.736 
^Mean of LG 1.76 and VG 1.75 

^Gited as Weiss, private communication. Using the GI value for the line strength from Weiss (1967) and the measured transition energy 
yields 1.77. 

^Weighted average of ten experimental values as of 2003. 


order many body perturbation theory [MBPT(4)] value of 71.7 calculated by Archibong and Thakkar |108] . while the 
coupled cluster with single and double excitation-effective Hamiltonian (CCSD-EH) approach of Stanton |117j yields 
72.2 using basis sets from Ref. Ell, denoted WMR in the Table. The pseudo-natural orbital coupled electron pair 
approximation (PNO-CEPA) calculation of Reinsch and Meyer |1()1) is close to the CI-I-MBPT calculation of Porsev 


and Derevianko m- Two time-dependent density functional theory (TDDFT) calculations are included in Table [HI] 


Using time-dependent DFT with a self-interaction correction (TDDFT-SIC), Chu and Dalgarno |109) obtained 71.8 
and using the symmetry adapted perturbation theory codes, SAPT(DFT), Patkowski et al. |119j obtained 73.27 for 
the polarizability. An extensive table of values for the polarizability of Mg resulting from various density functional 
theory functionals is given in Ref. E]. The R-matrix calculation of Robb |120j is 75(5); the relatively large value 
results because core-valence correlation effects were not included m- Excluding the relatively large value from 
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TABLE III: The static electric dipole polarizability (in units of Oq) for Mg in the ground state. Abbreviations for methods are 
defined in the text. 


Method 

«(0) 

Source 

Ref. (Year) 

Experiment (cryogenic molecular beam) 

59(15) 

Ma et al. 

[29] (2015) 

QDT 

69.54 

Chernov et al. 

[TTU] (2005) 

PP 

70.5 

Maeder and Kutzelnigg 

[90] (1979) 

CICP 

70.74(71) 

Muller, Flesch, & Meyer 

[TT3] (1984) 

CCD-bST (CCD) 

70.89 

Castro and Canute 

[TT^ (1993) 

Cl 

70.90 

Hamanou and Hibbert 

[52] (2008) 

MRCI 

71.2 

Partridge et al. 

[TTT) (1990) 

CICP 

71.35 

Mitroy and Bromley 

[97] (2003) 

CI-tMBPT 

71.3(7) 

Porsev and Derevianko 

[47] (2006) 

PNO-CEPA 

71.32 

Reinsch and Meyer 

PUTI (1976) 

MBPT(4) 

71.70 

Archibong and Thakkar 

[T08] (1991) 

TDDFT-SIC 

71.8 

Chu and Dalgarno 

[I09] (2004) 

MP 

72.0 

Patil 

\m\ (2000) 

MP 

72.1 

Victor and Slavsky 

[TTT| (1974) 

CCSD-EH (WMR) 

72.2 

Stanton 

[TTT] (1994) 

SAPT (DPT) 

73.27 

Patkowski et al. 

[TTn] (2007) 

CI-CV“ 

74.37 

Moccia and Spizzo 

[nS] (1988) 

Sum rule 

74.4(2.7) 

Reshetnikov et al. 

[HH] (2008) 

R-matrix* 

75(5) 

Robb 

[HO] (1975) 


“Mean of LG 74.7 and VG 74.03 

^Does not include core valence according to Ref. [90] 


Ref. [115) . the Cl, MRCI, CICP, and MBPT calculations fall in the range from 70.74 to 71.7. Reshetnikov et al. [89] 
use a semi-empirical method that utilizes a sum rule with constraints and error bars determined using measured 
the lifetime and excitation energies. Their value is 74.4(2.7), with the accuracy limited by the available input data. 
The recent experiment of Ma et al. |29j using a cryogenic molecular beam found a value of 59(15), which is not yet 
sufficiently accurate to test the calculations against. 

Table jlVj presents a collection of van der Waals constant values from the literature and a significant range is 
apparent, though the CI-I-MBPT and CICP calculations, which include models of core electron excitations, are in good 
agreement. The value 620(5) from the model potential (MP) calculation of Santra, Christ, and Greene |121j is close 
close to the pseudopotential (PP) calculation of 618.4 from Maeder and Kutzelnigg [90|, both of which include effective 
potentials to account for the presence of core electrons, but don’t fully include their excitations. The MP calculation 


from Patil [112j . however, is significantly larger, at 648. Three DPT calculations are listed in Table IV Hult et al. [122] 
introduced a local dielectric function and cutoff on the interaction volume and obtained 615. In contrast, Chu and 
Dalgarno used time-dependent DPT with a self-interaction correction (TDDPT-SIC) and an empirical correction to 
obtain 626 with an estimated uncertainty of 1%. Patkoswki et al. [119] used the symmetry adapted perturbation 
theory codes, SAPT(DPT), and obtained 635. The calculation of Cq by Stanton |117j used a quadrature and values of 
the dynamic polarizability at imaginary frequencies calculated using the CCSD-EH coupled cluster approach with the 
basis sets from |118| . The CI-CV calculations of Moccia and Spizzo [115] in the velocity gauge (VG) and in the length 
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TABLE IV: The dispersion constant ( 

are defined in the text. 

in units of 

e^flo) for the Mg dimer from various references. 

Abbreviations for methods 

Method 

Ce 

Source 

Ref. (Year) 

DFT 

615 

Hult at al. 

[Hg (1999) 

PP 

618.4 

Maeder and Kutzelnigg 

(HO] (1979) 

MP“ 

620(5) 

Santra, Christ, & Creene 

dH] (2004) 

TDDFT-SIC* 

626(6) 

Chu and Dalgarno 

[I09] (2004) 

CI+MBPT 

627(12) 

Porsev and Derevianko 

[13 (2006) 

CICP 

629.5 

Mitroy and Bromley 

[97] (2003) 

MP 

632.27 

Victor and Slavsky 

[TTT] (1974) 

SAPT(DFT) 

635 

Patkowski et al. 

[TT9| (2007) 

CCSD-EH (WMR) 

648 

Stanton 

[TTT] (1994) 

MP 

648 

Patil 

fm| (2000) 

CI-CV (VC) 

658.1 

Moccia and Spizzo 

fTT5| (1988) 

CI-CV (LC) 

670.9 

Moccia and Spizzo 

[145] (1988) 

Empirical 

683(35) 

Stwalley 

|HH] (1971) 

R-matrix 

689(70) 

Robb 

[Ho] (1975) 


“Does not include core contributions 

^Listed in Table VII as “corrected,” corresponding to an empirical rescaling. 


gauge (LG) are substantially larger than the other listed calculations. Robb |120j estimated his R-matrix calculation 
to be accurate to 10 %. Stwalley [SH] used an empirically constructed polarizability function to calculate Cq. The 
large value for Cg corresponds to the choice of 1.82 for the principal oscillator strength. In Ref. [5] it was found that 
a 2% uncertainty in Cg leads to an uncertainty of no more than 0.3 nm in the scattering length for ^^Mg 2 . 
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